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The  modulation  signal  for  the  rotating  helicopter  blade  is  derived 
in  Reference  1.  Two  basic  models  were  proposed  in  this  reference:  a 
model  based  on  scattering  from  the  blade  tip  and  a  specular  flash  model. 
Evidence  presented  in  Reference  2  is  somewhat  inconclusive,  although  the 
statement  is  made  that  "most  of  the  echo  is  derived  from  the  outer  portions 
of  the  rotor  blade  (approximately  the  outer  20%  of  the  blade  length)". 

Such  a  statement  supports  the  model  where  scattering  is  assumed  to  occur 
at  the  blade  tip. 

In  order  to  acccnmnodate  both  of  the  above  cases,  a  slightly  more 
general  model  for  the  modulation  signal  of  a  set  of  rotating  helicopter 
blades  than  given  in  References  3  and  4  is  defined  as 


V(t)  -V  A(0  )sinc[2(L  /X)sin0  ] 
X— #  n  ep  » 


J4»(L  /X)sin0 


where 


sinc(x)  -  (sinirx)/itx 


0n  -  2*<fst+n/N)  (2) 

fa  ■  spin  frequency 
N  •  number  of  blades 

L _  ■  Lcosa,  L_  -  effective  blade  length 

•p  #  c 

L _  ■  L.coaa,  L_  -  location  of  phase  center 

cp  c  c 

a  ■  angle  between  LOS  and  plane  of  rotation 
X  ■  wavelength 
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The  quantity  A(6  )  is  used  to  allow  different  scattering  properties  on 
the  leading  and  lagging  edges  of  the  blade.  Since  0q  is  referenced  to 
the  time  of  occurrence  of  the  maximum  Doppler,  we  can  define 

{1  for  leading  edge  (cos8^>0) 

(3) 

p  for  lagging  edge  (cos0q<O) 

The  Instantaneous  phase  and  frequency  of  the  phase  center  of  the  n**1  blade 
are  given  by 

♦i  -  4s<Lcp/X)sin0n  (4) 

f±  -  4w(Lcp/X)fscos0n  (5) 

where  0  is  given  by  (2). 

Q 

From  (5)  we  note  that  the  maximum  Doppler  is  given  for  cosS^  *  1  as 

f  "  4w<LPn/X>fa  <6> 

max  cp  s 

In  a  typical  case  we  might  have  Lcp  “4  m,  X  “  .02  m,  and  fg  ■  10  Hz,  so 

that  f _ “  25.1  kHz.  One  half  of  a  cycle  later  the  minimum  Doppler  would 

be  -f^,  so  the  bandwidth  of  the  modulation  signal  is  2fma^,  or  50.2  kHz 
for  the  example.  In  order  to  use  (1)  to  simulate  the  modulation  signal, 
the  sampling  frequency  must  be  at  least  as  large  as  2f  to  prevent 
aliasing  or  foldover  (the  amplitude  modulation  in  the  sine  function  in 
(1)  causes  some  additional  Doppler  broadening  beyond  ±ftway) ♦  However,  in 
the  proposed  application  this  would  stipulate  a  sampling  frequency  of  at 
least  50  kHz,  which  is  about  a  factor  of  4  larger  than  the  receiver  pro¬ 
cessing  bandwidth.  In  other  words,  (1)  simulates  signals  that  would  be 
outside  of  the  receiver  band,  as  well  as  the  desired  in-band  signal. 
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In  order  to  minimize  the  computation  load  it  would  be  desirable  to 
sample  at  the  rate  that  is  coincident  with  the  receiver  processing  band. 
But  in  order  to  prevent  the  folding  over  of  frequencies  outside  of  the 
desired  band,  we  must  suppress  the  unwanted  signals.  This  will  be 
accomplished  by  computing  (5),  the  instantaneous  frequency  of  each 
component,  to  determine  if  that  component  is  within  the  receiver  pro¬ 
cessing  band.  If  it  is,  then  that  component  will  be  included  in  the 
summation  in  (1) . 

Disregarding  the  return  from  the  hub  portion  of  the  blade  assembly 

for  the  moment,  it  can  be  shown  that  as  long  as  the  number  of  blades  is 

relatively  small  (s  6),  the  maximum  number  of  components  of  (1)  that  can 

fall  within  a  band  that  is  251  of  2f  is  two.  Thus  for  N-2  of  the 

max 

components,  only  f ^  in  (5)  need  be  computed.  This  result  significantly 

reduces  the  computation  time  for  N  >  2. 

The  hub  portion  of  the  blade  assembly  also  scatters  energy.  We  can 

simulate  this  component  of  the  modulation  signal  by  the  use  of  (1)  if  we 

redefine  L  and  L  .  In  general,  the  scattering  will  be  much  more  iso- 
ep  cp 

tropic  than  for  the  blade  so  that  L  will  be  small  (if  L  -  0  then  it 

ep  ep 

would  be  isotropic).  If  the  spectral  band  for  the  hub  signal  falls  within 
the  desired  processing  band,  then  occasionally  the  blade  Doppler  will  also 
fall  within  this  band.  But  this  occurs  only  when  the  blade  is  near  the 
quadrature  points  to  the  maximum  and  minimum  Doppler.  In  other  words, 
we  will  be  looking  at  the  ends  of  the  blade.  In  order  to  conserve  computer 
resources,  the  assumption  is  made  that  when  the  hub  return  falls  within 
the  processing  band,  the  blade  return  will  be  neglected. 

References 

1.  Mitchell,  R.  L.,  and  I.  P.  Bottllk,  "Techniques  for  Simulating 
Realistic  RF  Environment  Signals  on  the  RFSS,"  MRI  Report  131-25, 

28  February  1977. 

2.  "Improved  Hawk  Systems  Threat  Definition  and  Capabilities  Against 
Helicopter  and  Liaison  Aircraft,"  Raytheon  Final  Report  BR-6501, 
Revision  B,  15  Jan  1975. 

3.  Mitchell,  R.  L.,  "Helicopter  Blade  Modulation  Model  (Revised)," 

MRI  Report  149-16,  29  March  1978. 

4.  McPherson,  D.  A. ,  "Rotating  Blade  Modulation  Waveform,"  Boeing/ 
Huntsville  Report,  15  March  1978. 


4 


\ 

Appendix 
Program  Listings 


Routine  Comments 

TEST  Test  program 

HBMSD  Set-up  parameters  for  call  to  HEMS. 

Probably  will  be  done  in  Datacraft/6. 

HEMS  Generates  modulation  signal. 

Probably  will  be  done  in  AP120B. 

Fortran  version  of  subroutine  to  be 
implemented  as  table  lookup. 


SINCOS 


PROGRAM  TEST (OUTPUT) 

COMMON  /Tl/  VR  <  200 )  x 

COMMON  /T2/  VI (200) 

COMMON  /Dl/  NSAMP,  FR,  FS,  FO,  Fl»  R,  DELH,  WL 

COMMON  /D2/  XFS,  XF1,  XF2,  ASK,  ABL,  NBL.  XLEP,  XLCP,  RHO.  PEFF 

NSAMP-200 

FR- 10000. 

FS“10. 

FO-O. 

FI— 5000. 

R-1000. 

DELH-O. 

WL«.  02 
CALL  HBMSD 

PRINT  100,  <K, VR(K), VI <K>, K=l,  NSAMP) 

PRINT  101,  N3L.  XFS,  XF1,  XF2.  ASK,  ABL,  XLEP,  XLCP.  RHO,  PEFF 
100  FORMAT  ( 1 6,  2F 1 2.  6 ) 


101  FORMAT </I6,  SF12.  6,  E16.  6) 


SUBROUTINE  HBMSD 


C  ' 

C  DRIVER  FOR  HBMS 
C 

C  IN  THIS  SUBROUTINE  WE  SET  UP  THE  PARAMETERS  FOR  ONE  CALL  TO  SUBROUTINE 
C  HBMS,  WHICH  GENERATES  THE  HELICOPTER  BLADE  MODULATION  SIGNAL. 

C 

C  THE  REAL-TIME  INPUTS  THROUGH  COMMON  /Dl/  ARE . 

C 

C  NSAMP  «  NUMBER  OF  SAMPLES  COMPUTED  .PER  CALL  TO  HBMS 
C  FR  -  SAMPLING  FREQUENCY 

C  FS  *  BLADE  SPIN  FREQUENCY 

C  FO  *  DOPPLER  FREQUENCY  OF  SKIN  RETURN 

C  FI  *  LOWER  FREQUENCY  OF  RECEIVER  PROCESSING  WINDOW 

C  R  “  RANGE 

C  DELH  -  DIFFERENCE  IN  ALTITUDE  BETWEEN  RADAR  AND  HELICOPTER 
C  WL  -  WAVELENGTH 

C 

C  THE  OUTPUTS  THROUGH  COMMON  /D2/  ARE . 

C 

C  XFS  -  FS/FR 

C  XF1  ■  (FI— FO) /FR 

C  XF2  -  XF1+1. 

C  ASK  ■  REFERENCE  AMPLITUDE  OF  SKIN  RETURN 

C  ABL  -  REFERENCE  AMPLITUDE  OF  BLADE  OR  HUB  RETURN 

C  NBL  -  NUMBER  OF  BLADES 

C  XLEP  -  NORMALIZED  PROJECTED  EFFECTIVE  BLADE  LENGTH 

C  XLCP  ■  NORMALIZED  PROJECTED  EFFECTIVE  LOCATION  OF  PHASE  CENTER 

C  RHO  ■  RATIO  OF  BLADE  AMPLITUDE  ON  LEADING  EDGE  TO  LAGGING  EDGE 

C  PE FF  -  EFFECTIVE  POWER  TO  BE  RADIATED  FROM  RFSS  ARRAY 

C 

C  ALL  OF  THE  REMAINING  PARAMETERS  ARE  NON-REAL -TIME  AND  ARE  SET  IN  THIS 

C  SUBROUTINE.  FIRST  WE  HAVE . 

C 

C  AS  *  AMPLITUDE  OF  SKIN  RETURN  (=SQRT(RCS>,  NONFLUCTUATING) 

C  PTODSQ  ■  PRODUCT  OF  TRANSMIT  POWER,  GAIN,  SQUARE  OF  CHAMBER  LENGTH 
C 

C  THE  BLADE  PARAMETERS  ARE . 

C 

C  NBLB  -  NUMBER  OF  BLADES 

C  ALEB  -  EFFECTIVE  BLADE  LENGTH 

C  ALCB  -  EFFECTIVE  LOCATION  OF  PHASE  CENTER 

C  RHCB  -  RATIO  OF  LEADING  TO  LAGGING  EDGE  AMPLITUDE  FOR  BLADE 

C  AB  -  PEAK  AMPLITUDE  OF  BLADE  RETURN  (=SQRT(RCS>) 

C 

C  THE  HUB  PARAMETERS  ARE . 

C 

C  NBUH  -  NUMBER  OF  HUB  SECTIONS  (USUALLY  EQUAL  TO  NBLB) 

C  ALEH  ■  EFFECTIVE  BLADE  LENGTH  FOR  HUB  (PROBABLY  SMALL) 

C  ALCH  *  EFFECTIVE  LOCATION  OF  PHASE  CENTER  FOR  HUB 


ooo  ooo  ooo  ooo  ooooooo 
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RHOH  *  RATIO  OF  LEADING  TO  LAGGING  EDGE  AMPLITUDE  FOR  HUB 
AH  -  PEAK  AMPLITUDE  OF  HUB  RETURN  <=SGRT<RCS>) 

THE  TARGET  IS  ASSUMED  TO  BE  NONFLUCTUATING. 

THE  UNIT  OF  DISTANCE  MUST  3E  CONSISTENT  THROUGHOUT. 

COMMON  /Dl/  NSAMP*  FR.  FS»  FO*  FI*  R»  DELH.  WL 

COMMON  /D2/  XFSi  XF1.  XF2.  ASK.  ABL.  NBL.  XLEP.  XLCP.  RHQ.  PEFF 

DATA  NBLHi  ALEH,  ALCH.  RHOH.  AH/3.  0.  .  .  5.  1.  ,  1.  / 

DATA  NBLB.  ALEB.  ALCB.  RHOB.  AB/3.  2.  .  8.  .  .  3.  1.  / 

DATA  AS.  PTGDSG/1.  .  1.  / 

DATA  FOURPI/12.  566370614/ 

\ 

COMPUTE  COSINE  OF  ASPECT  ANGLE 

CA-SQRT <  1 .  -  < DELH/R ) **2 ) 

COMPUTE  NORMALIZED  FREQUENCIES  . 

XFS-FS/FR 

XFO-FO/FR 

XFl-Fl/FR-XFO 

XF2-XF1+1. 

DETERMINE  IF  HUB  DOPPLER  IS  VISIBLE 

IF(XF1.  Lt.  0.  .  AND.  XF2.  CT.  0.  )  GO  TO  20 

COMPUTE  BLADE  PARAMETERS 

NBL-NBLB 
XLEP«ALEB*CA/WL 
XLCP«ALCB*CA/WL 
RHO-RHOB 
APEAK-AS+2.  *AB 
ASK-AS/APEAK 
ABL-AB /APEAK 
00  TO  30 

COMPUTE  HUB  PARAMETERS 

20  NBL-NBLH 

XLEP-ALEH*CA/WL 
XLCP-ALCH*CA/WL 
RHO-RHOH 
APEAK-AS+2.  «AH 
ASK-AS/APEAK 
ABL-AH/ APEAK 


30  PEFF=PTGDSQ* ( APEAK**2 ) / <  FOURP I #R#*4 
CALL  HBMS  ' 

RETURN 

END 


ooooooooooooonooooaooooooooooooooo 


SUBROUTINE  HBMS 

X 

HELICOPTER  BLADE  MODULATION  SIGNAL 

IN  THIS  SUBROUTINE  WE  COMPUTE  NSAMP  COMPLEX  SAMPLES  IN  THE  ARRAY-PAIR 
(  VR»  VI )  OF  THE  HELICOPTER  BLADE  MODULATION  SIGNAL  AT  THE  SAMPLE  RATE 
WHICH  WILL  BE  DESIGNATED  AS  FR.  ALL  INPUT  FREQUENCIES  ARE  NORMALIZED 
TO  THIS  QUANTITY,  NAMELY 

'  XFS  *  FS/FR  (FS  =  BLADE  SPIN  FREQUENCY) 

XF1  *  Fl/FR  (FI  =  LOWER  FREQUENCY  OF  PROCESSING  WINDOW) 

XF2  *  F2/FR  (F2  =  UPPER  FREQUENCY  OF  PROCESSING  WINDOW) 

FI  AND  F2  ARE  REFERENCED  TO  THE  SKIN  DOPPLER  (DC). 

WE  ASSUME  THAT  THE  PROCESS  IS  UNDERSAMPLED  SO  THAT  ONLY  THOSE  DOPPLER 
FREQUENCIES  BETWEEN  FI  AND  F2  WILL  BE  SIMULATED  (IF  ALL  FREQUENCIES 
WERE  TO  BE  SIMULATED  THEN  FOLDOVER  WOULD  RESULT  WHENEVER  8*PI*XLCP*FS 
EXCEEDED  F2-F1).  IN  ADDITION,  WE  HAVE  ON  INPUT 

XLEP  =  LEP/WL  ( LEP=EFFECT I VE  PROJECTED  BLADE  LENGTH) 

XLCP  =  LCP/WL  (LCP=EFFECTIVE  PROJECTED  LOCATION  OF  PHASE  CENTER) 
ASK  ■  REFERENCE  AMPLITUDE  FOR-  SKIN  RETURN 
ABL  »  REFERENCE  AMPLITUDE  FOR  BLADE  RETURN 

RHO  *  RATIO  OF  BLADE  AMPLITUDE  ON  LEADING  EDGE  TO  LAGGING  EDGE 
NBL  -  NUMBER  OF  BLADES 

NSAMP  =  NUMBER  OF  SAMPLES  SIMULATED  PER  CALL  TO  HBMS 
THE  BLADE  LENGTH  REFERS  TO  THE  RADIUS  OF  THE  ORBIT. 

IN  THE  USE  OF  THIS  SUBROUTINE  IT  IS  ASSUMED  THAT  XLEP  AND  XLCP  WILL 
DESIGNATE  THE  PARAMETERS  OF  THE  HUB  PORTION  OF  THE  BLADE  ASSEMBLY  IF 
FI  AND  F2  ENCOMPASS  DC  (OR  THE  DOPPLER  OF  THE  SKIN  RETURN).  OTHERWISE 
THE  PARAMETERS  WILL  BE  FOR  THE  ACTUAL  BLADE. 

COMMON  /Dl/  NSAMP 

COMMON  /D2/  XFS,  XF1,  XF2,  ASK,  ABL,  NBL,  XLEP,  XLCP,  RHO,  PEFF 
COMMON  /Tl/  VR(1) 

COMMON  /T2/  VI ( 1 ) 

DATA  TWOPI/6.  283185307/ 

DATA  PHS/O.  / 

XN-1.  /FLOAT (NBL) 

XLCP2-XLCP*2. 

X 1 -TWOP I *XLCP2*XFS 

X2-TW0PI*XLEP 

DO  30  1*1,  NSAMP 

UR-ASK 

UI-O. 

DO  20  L*l, NBL 

CALL  8INC0S(S1,  Cl,  PHS) 


L 


XFI»X1*C1 

IF(XFI.  LT.  XF1 )  GO  TO  15 
IF(XFI.  GT.  XF2)  GO  TO  15 
AR0-XLCP2*S1 
CALL  SINC0S<S2,  C2,  ARG) 

ARC-X2*S1 

A-ABL 

1F( ABS< ARG).  GJ.  .  001)  A=ABL*SIN< ARG)/ARG 
IF(C1.  LT.  0.  )  A=RHO*A 
UR«UR+A*C2 
UI«UI+A*S2 
15  PHS=PHS+XN 
80  CONTINUE 
VR<I)=UR 
VI<I)=UI 
PHS-PHS+XFS 
30  CONTINUE 

PHS*PHS— FLOAT ( NSAMP ) 

RETURN 

END 


SUBROUTINE  SINCOSCS, C. X) 
DATA  T WOP  1/6.  2B3185307/ 
S»SIN<TWOPI*X> 
C-COS<TWOPI*X> 

RETURN 
END 


